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Abstract 

We study the nonlinear hopping transport in one-dimensional rings and open channels. Analyti- 
cal results are derived for the stationary current response to a constant bias without assuming any 
specific coupling to the external fields. It is shown that anomalous large effective jump lengths, 
as observed in recent experiments by taking the ratio of the third order nonlinear and the linear 
conductivity, can occur already in ordered systems. Rectification effects due to site energy dis- 
order in ring systems are expected to become irrelevant for large system sizes. In open channels 
in contrast, rectification effects occur already for disorder in the jump barriers and do not vanish 
in the thermodynamic limit. Numerical solutions for a sinusoidal bias show that the ring system 
provides a good description for the transport behavior in the open channel for intermediate and 
high frequencies. For low frequencies temporal variations in the mean particle number have to be 
taken into account in the open channel, which cannot be captured in the more simple ring model. 

PACS numbers: 66.30.H-,05.60.Cd,66.10.Ed 
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I. INTRODUCTION 



The particle transport in one-dimensional systems is of vital interest for many problems in 
physics and biology A prominent example is the electron or hole transport in the operation 
of conducting nanowires, including molecular wires.- In such systems transport can be dom- 
inated by quantum mechanical tunneling or band motion (the coherent transport limit) but 
many systems belong to the hopping transport limit, where conduction is a manifestation of 
succession of many incoherent hopping steps.- 1 ^ For example, both conduction mechanisms 
were observed in different DNA sequences.- One-dimensional hopping motion is also the 
decisive transport mechanism in ion conduction through membrane channel** 1 ^ and unidi- 
rectional motion of motor proteins along filaments. 8 ' 9 In the connection of the latter example 
much attention have received recently boundary driven phase transitions in one-dimensional 
lattice gases with site exclusion and asymmetric hopping dynamics, commonly referred to 
as "asymmetric site exclusion process" (ASEP), or, in case of unidirectional transport, as 
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14. Re- 



"totally asymmetric site exclusion process " (TASEP) - for reviews, see Refs. 
cently, properly modified models^ were applied to describe the transport of single-stranded 
DNA segments through nano channels.™ 

The treatment of one- dimensional systems is moreover frequently used as a starting point 
for describing transport processes in higher dimensions, since it often allows one to derive 
analytical results. In transferring essential results to higher dimensions one has, however, 
to be careful. An example is the tracer diffusion in one-dimensional hard-core lattice gases, 
which exhibits a subdiffusive behavior for long times that originates from the fact that 
particles cannot pass each other in one dimension.— 

In this work we will study the thermally activated hopping conduction in one-dimensional 
lattices for non-interacting particles in arbitrary energy landscapes. In particular we consider 
the nonlinear transport in strong static and periodic fields. For couplings oc exp(±w/2) 
of the external bias u to the bare hopping rate, this problem was first studied for ring 
systems (periodic boundary conditions) in Ref. isl An exact result for the stationary current 
was derived, in generalization of an analogous treatment for Brownian dynamics. 19 As a 
particularly interesting feature, rectification effects were shown to be present for energy 
landscapes with site energy disorder. 

The problem got renewed interest recently for describing measurements on thin glassy 
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electrolytes under high voltages,— i2ii22i22 which allow one to reach the weak nonlinear regime, 
u = qEa/k-gT ~ 1, where q is the charge of the mobile ions, a is a typical hopping distance 
of 2-3A, E is the applied electric field, and k^T the thermal energy. In these experiments no 
rectification were observed so far, meaning that the current turned out to be an odd function 
of the applied field. On the other hand, these measurements can be used to determine 
an effective length scale a e fr when analyzing the ratio 03/01 of the third order nonlinear 
conductivity 03 to the linear conductivity <j\ (cf. Eq. [T7] below). This length a c s appears 
to be unphysically large if it is compared to typical jump lengths a ~ 2 — 3 A. Such 
comparison is motivated by the result^ j dc oc smh(qEa/2kBT), which applies to the most 
simple situation of single-particle hopping in an ordered system with the aforementioned 
coupling oc exp(±w/2) of the bias to the bare hopping rates (see below). For different glassy 
electrolytes a e g either increases or decreases with T (in the temperature ranges studied a 
linear behavior was observed). It was also found that 03 > 0, while a 5 has different sign 
for different glass compositions. In the frequency-dependent response the real part j' 3 {u>) 
of the third order harmonics jz{uj) has a negative sign for low frequency. With increasing 
frequency, j 3 (^) increases and becomes positive close to the onset frequency of the dispersive 
part in the first order harmonics j[{uj) (which gives the linear response conductivity a\ (a;)). 
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Taking disorder averages 2 - of the analytical expression for the current derived in Ref. 
was suggested that the large values of a e g have their origin in the spatial variation of hopping 
rates in the glassy material. Moreover, based on a small u expansion, it was predicted that 
a e fr oc iV 1//2 , where iV ~ L/a is the number of sites of the film sample in field direction. 
However, this result followed when expanding terms as exp(Nu) in the analytical result for 
the current in power of Nu. Since the nonlinear transport becomes relevant for u ^ 1, 
and N should be significantly larger than one (to avoid boundary effects), this expansion 
in powers of Nu is in general not appropriate. Rather one should take the thermodynamic 
limit N — > 00 before carrying out the small u expansion of the current ; 25 - which can yield 
non-analyticities in the current response. It was argued 2 ^ that these non-analyticities could 
spoil the analysis of nonlinear conductivities based on odd powers in the field amplitude, as 
they are commonly employed in experiments. 

An open question is whether the rectification effects occurring in finite systems are present 
also in the thermodynamic limit. Intuitively, one would expect that in the absence of long- 
range correlations in the energy landscape (i.e. correlations decaying faster than 1/distance), 
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self- averaging effects suppress rectification properties the more the larger the system size 
becomes. As a consequence one would predict rectification effects to disappear in the ther- 
modynamic limit. While this is in agreement with experimental observations (for sample 
thicknesses so far studied), it has not yet been demonstrated by theoretical analysis. To 
avoid the problem of possible rectification effects and to enforce that the current is an 
odd function of u, energy landscapes with point symmetry were considered in Refs. 
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However, the constraint of point symmetry implicitly introduces long-range correlations in 
the energy landscape and it is questionable if such procedure is suitable to describe real 
experimental situations. 

In this work we will treat the following open problems: 

(1) Analytical results for the stationary current in ring system with M sites were derived 
up to now for the coupling cx exp(±ii/2) of the external bias u = qEa/ksT to the bare 
rates (rates in the absence of the external driving). This rate emerges naturally when 
approaching the hopping limit of the overdamped Brownian dynamics (Smoluchowski 
equation) of noninteracting particles. However, in interacting many-particle systems 
more complicated couplings of the rates to the external field can be imagined, when 
mapping the dynamics to an effective one-particle hopping process in a renormalized 
energy landscape. We therefore derive the stationary current for arbitrary couplings, 
and discuss in more detail the behavior for jump rates obeying the condition of detailed 
balance. We find that is then possible to obtain already in an ordered system effective 
lengths scales a e ^ significantly larger than the jump length a. Hence it appears that not 
only the disorder affects a e ff. 

(2) As outlined above, for relating the theoretical results to experiments in the nonlinear 
regime, one should first perform the thermodynamic limit M — > oo before expanding the 
current in powers of the field amplitude. By performing this limit we also clarify the role 
of rectification effects for large M. 

(3) For ring systems it is unclear how the periodic boundary conditions affect the stationary 
current. We therefore study the analogous problem in an open channel, where particles 
are injected and ejected from two particle reservoirs on the left and right side with 
electrochemical potentials /xl and //r, respectively. The rates for the local exchange of 
particles with the reservoirs fulfill detailed balance with respect to the grand-canonical 



ensembles associated with «l and /jr. We will treat the linear limit of the rate equations 
in this work to avoid boundary induced phase transitions as occurring in ASEPs or 
TASEPs.^^ 

(4) Up to now the time-dependent nonlinear current response has been rarely studied. 39 
Here we will investigate for both the ring systems and open channels this time-dependent 
response to a sinusoidal driving with large field amplitude Eq by numerically solving the 
corresponding rate equations for the occupation probabilities. The data are analyzed, by 
using standard Fourier analysis, in terms of harmonics j n (uj) of nth order. We present 
results for spatially uncorrelated barrier energies with uniform distributions and discuss 
the relation of the harmonics in the ring and open channel with respect to different 
frequency regimes. 



II. TRANSITION RATES AND ENERGETIC DISORDER 

For convenient notation, we define k-^T as the energy unit in the following, k^T = 1. 
In a disordered energy landscape with site energies and energy barriers £4,fc+i = £4+i,fc 
between sites k and k + 1, the rates T^(t) and T^ +1 (t) are considered to be functions of e^, 
€k+i, and U^k+i- In addition they depend on the external bias u(t), which we assume to 
be homogenous over the ring or channel, corresponding to a linear decrease of the external 
potential. If the rates obey detailed balance at each time instant, their ratio rjk{t) is given 
by 

r+(t) 

Vk(t) = = exp(-AE M+1 ) , (1) 

where 

A-E fcjfc+ i = e k+1 -e k -u. (2) 

In the presence of screening effects, the assumption of a constant potential gradient is not 
valid, leading to a bias depending on k. The analytical formulae derived in the following 
sections can be generalized to this situation. 

To illustrate our findings we will consider two types of rates and two types of energetic 
disorder. For the rates these are the "exponential rates" 

= -exp(-C/ fejfc+ i)exp(-A£; fejfc+ i/2) (3) 
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and the Glauber rates 26 

r t = |exp(-C/ fc)fe+ i) 

where 7 is a bare jump rate. For the energetic disorder, we consider either pure barrier 
disorder (all e k = 0), or pure site energy disorder (all barriers U kjk+ i — 0). The barrier 
and site energies are uncorrelated random variables drawn from box distributions, U k>k+ \ £ 
[0, Ajj] and 6i G [— A e /2, A e /2] with widths Ajj and A e , respectively. 

III. CONDUCTION IN RING SYSTEMS 

We study the nearest neighbor hopping of one particle on a ring with M sites % — 1, . . . , M. 
The rates for a jump from site i backward and forward at time t are denoted as T~(t) and 
r^"(t), respectively. The probabilities Pi(t) for the particle to be on site i at time t obey the 
rate equations 

Pi =jir-i,i(t) -j i>i+1 (t) , i = l,...,M, (5) 

with the local currents 

ji, i+1 (t) = r+(t) Pi (t) - rr +1 (t) Pi+1 (t) , i = 1, . . . , M . (6) 

In writing Eqs. ( |5|6l) and further equations below we implicitly assume that the periodic 
boundary conditions are taken into if the index i falls out of the range 1, . . . , M, i.e. pi + M{t) = 
Pi{t), Ff +M (t) = Tf(t), j i+M ,i+i+A4(t) = etc. The rate equations preserve the 

normalization YliLiPi(t) = 

Due to the normalization of the occupation probabilities to one particle, the current j 
refers to the single particle current. If we consider a fixed number density n per lattice site 
of non-interacting particles, the total current is 

J = nMj . (7) 

In the case of charged particles the corresponding charge current per lattice site is qJ and 
the charge current density qJ/A, where A is a cross sectional area associated with each 
lattice bond. 



1 + tanh 



/ AE kk < 
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exp(-C4 |fe+ i) 
1 + exp(-AE kk+ i) 



(4) 
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A. DC current 



In a static (time-independent) driving field u the system reaches a stationary state for 
long times, where the occupation probabilities become constant, p^ = pf, and all local 
currents in Eq. (j^D are equal, j^+i = jdc- Setting = 1/T^ +1 this leads to the recursion 
relation 

pT+i = ViPf - Ktfdc, (8) 

with solution 

EM T-r/c :-l 
, k=l K i-k 11;=1 Vi-l ( Q \ 

Pi — Jdc M ■ \?) 

U k =i Vi-k - 1 

The current jdc follows from the normalization, 

j- = Z n? ?~ l > (10) 
l= i U k =i Vi-k - 1 

which in turn fixes the occupation probabilities pf. Equations (MTUj) hold true for arbitrary 
set of rates (as long as they do not exclude the formation of a unique stationary state). 

For detailed balanced rates these expressions can be simplified. With condition we 
have 

k 

| J Vi-l = exp(ei_ fc - + ku) (11) 
i=i 

so that Eq. ffTUl) can be written in the form 

- - 4^7 1 <f t ' + Ct+1>/2 ' «P(~^) ■ M 



3dc e ~ 1 *=i Vr+(«)r- +1 («) 

where we explicitly indicated the dependence of the jump rates = T^(-u) on the external 
bias u. For the coupling oc exp(±ii/2) of the rates to the external field it can be shown that 



this formula agrees with Eq. (10) in Ref. Il8| (or with Eqs. (8-11) in Ref. l22l ) . 

In the linear response limit u — > 0, Eq. (112"]) reduces to the result^ qJ& c /A = qnMj/A = 
<J\Eq with 

ffI =fer(Mg^J ■ (13) 

where p e ^ cx exp(— e&) is the equilibrium distribution and are the rates in the absence 
of external driving (u = 0). This formula can be viewed as resulting from conductances 
oc p e ^T^ in serial order. 



In systems with only barrier disorder (all sites have the same energy = 0), Eq. (|T2|) 
reduces to 



Because r^(— u) = r fe+1 (w) in this case, we have T^(— n)r fc+1 (n) = r fc+1 (M)r^" (w) = 
T~l(u)r^(u), and it follows that j(—u) = —j{u) for each disorder configuration. This is 
at first sight a surprising results, since one could consider an asymmetric spatial arrange- 
ment of barriers, for example, U^i+i = iUo for i = 1, . . . , iV and Uo > 0. If a particle would 
be driven in the direction of increasing i, it encountered increasing barriers until a jump from 
the largest to the smallest barrier occurs (after passing the barrier Un,n+i = E/jy,i between 
sites N and 1). When driving the particle in the reverse direction the opposite behavior 
would results, i.e. the particle encountered smaller and smaller barriers until a jump from 
the smallest barrier to the largest occurs. 

Moreover, as long as the barriers for the local transitions are taken into account by a 
simple Boltzmann factor, i.e. T^(m) oc exp(— Uk,k+i)f+(u) and Y^ +l {u) oc exp(— Uk,k+i)f-{u) 
with functions f±(u) independent of k, one obtains the same current- voltage curve J{u) = 
nMj(u) as in an ordered system up to a rescaling factor. In such ordered system, T^r^ +1 = 
T + r _ is independent of k, and one obtains an M independent total current J& c = nMj dc , 



Figure [T] shows the current Jdc(w) in the ordered ring system (or in the ring systems with 
barrier disorder) for the exponential rates ([3]) and the Glauber rates (j3J). For comparison we 
also show the average current J<i c (u) in the case of the box distribution of site energies with 
A e = 6. The current was calculated according to Eqs. (ITfTiZl) and averaged over 10 3 different 
realizations of the site energy disorder in rings with M = 10 3 sites. The current-voltage 
curves in the experimentally relevant regime u ^ 1 tend to have a more convex shape in the 
presence of site energy disorder. For the Glauber rates the current is smaller and saturates 
for u — * ±oo. 

In the ordered ring system, for the generic coupling r ± (u) = (7/2) exp(±w/2), one re- 
covers from Eq. (fl5il the known result for the charge current density 24 
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FIG. 1: Current jd c (u) in the ring system as a function of the bias u for the exponential rates 
(solid line) and the Glauber rates (dashed line). Results are shown for an ordered system and a 
box distribution of site energies with A e = 6. In the system with site energy disorder the mean 
current is shown, obtained after averaging jd c from Eq. (llOj) over 100 realizations. In the case 
of pure barrier disorder the same curves as in the ordered system are obtained for each disorder 
realization up to a (realization-dependent) rescaling of the current (see text). 



with (Ji = (n/ 'aA) r yq 2 a 2 jlk^T and 03 = [nj l aA) r yq 4 'a 4 ' / '^(hsT) 3 These results motivate to 
define an effective jump length by 



24(AfeT) 2 <7 3 



(17) 



<r 0-1 

However, even in an ordered system it is possible that this effective jump length does not 
yield a reasonable estimate of the true jump length a. The reason is that, while the linear re- 
sponse quantity o\ is universal (i.e. independent of the specific form of the jump rates), this is 
not the case for the nonlinear conductivity 03. For example, for the Glauber rates, we obtain 
J = 7ntanh('u/2) from Eq. (j!5p and accordingly a negative 03 = — (n/aA)'yq 4: a 4: /24(kBT) 3 . 
If this would be inserted in Eq. (|T7|) . a e s became imaginary. 

In the general case, we can expand T + (u) in a Taylor series, T + (u) = ( r y/2)(l + a\u + 
ct2ii 2 + a^u 3 + ...). From T~(u) = T + (—u) = exp(— u)T + (u) it follows that a\ = 1/2 
independent of the specific form. With Eq. ( |T5|) we find 0-3/0-1 = («2 — l/12)g 2 a 2 /(fcB^) 2 , 
i.e. 



l eft 



24 ( « 2 - 1 ] a' . 



(18) 



We conclude that dependent on a 2 (e.g., «2 = 1/8 for the exponential rates, yielding a e s = a, 
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and a 2 = for the Glauber rates, yielding a 
in an ordered system. 



cff 



-2a 2 ) different a e g can be obtained even 



B. Thermodynamic limit and rectification 

In the thermodynamic limit M — > oo the sum over k in Eq. (T12I) can be replaced by 
a disorder average (...) if the site energies and energy barrier do not exhibit very broad 
distributions or long-range correlations, i.e. if the the system is self-averaging. Accordingly 
we define 

M 



aAu) 



lim 



exp[(e fc + e fc+i )/2] 



k=l 



CXp(-€ k+l ) 



exp[( ei - e 2 )/2] 



y / T + (u; U 12 , ei, e 2 )T-(u] U 12 , ei, e 2 ) 

exp[(ei + e 2 )/2]exp(-e 3 ) 
\/T+(u; U 12 , ei, e 2 )r-(w; C/ 12 , e l} e 2 ) 



M-l 



(19) 



exp[(e 2 - ei )/2] 



Z = M . 



T+( U ;C/ 12 ,e 1) e 2 )r-( U ;r7 12) e 1) e 2 ) 

where we took into account the periodic boundary conditions and have explicitly denoted 
the dependence of the jump rates T^(u) = r ± (w; Uk,k+i, ^k, e fe+i) on the energies. 

Keeping the number density n fixed in the limit M — ► oo, we then obtain from Eq. (|12p 
for the total current 

>l 



^dc(^) = 2n sinh 



(20) 



• 2/ [0(-ii) ai («) + ^(n)a M ( M )](e |M| - 1) + a 2 (u) ' 

where 0(.) is the Heaviside step function [6(x) = 1 for x > and zero else]. As discussed 
in the Introduction, Eq. (120]) should apply to typical experiments on thin film electrolytes. 
There should be no notable dependence of the current (and the nonlinear conductivities) on 
the film thickness, in agreement with the experimental observations. 

We can further show that the current from Eq. ( 1201) is anti-symmetric with respect to the 
bias u. To this end we have to analyze the symmetry properties of the ai(u). Note that in the 
averages in Eq. (Tl9~|) there occur configurations with two or three sites only, having mutually 
independent random site energies ei, . . . , e^. As illustrated in Fig. [2], to each realization of the 
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two energies e\ and e 2 there exists a "mirror configuration" with interchanged site energies e\ 
and e 2 , and the same value of 63. Since these mirror configurations occur with equal statistical 
weight and exhibit the symmetry property T + (u; U% 2 , ei, e 2 ) = r~(— u; U% 2 , e 2 , e%), we can 
use T + (u;Ui2,ei,e 2 )T'(u;Ui2,e 1 ,e 2 ) = T + (-u; U X2 , e 2 , e 1 )Y~{-u] U 12 , £2, <n) in the averages 
of Eq. CHI . This implies ai(— w) = a^{u) and a 2 (— w) = a 2 (w), leading to J(—u) = —J(u). 
Let us note that this does not imply that the expansion of J(u) contains odd powers of u 
only. Terms cx |ii| 2n+1 w, n — 0, 1, . . ., can occur according to Eq. (|20|) (see also the discussion 



in Ref. 



23l for the consequences of these non-analytic terms with respect to the analysis of 



experiments). 



-u 





FIG. 2: Two mirror configurations with interchanged site energies e\ and e 2 (and same €3, not 
shown), as appearing with equal statistical weight in the averages in Eq. (|19p . The jump rate 
T + (n, U\2, €±, e 2 ) in the left configuration is equal to the jump rate T~(—u, U12, e 2 , ei) in the right 
configuration after reversal of the bias u. 



In view of the antisymmetric current in the thermodynamic limit, we expect, due to 
self- averaging, rectification effects for one system to become smaller with increasing system 
size. To check this expectation, we define the rectification parameter 

J dc (u) + Jdc(-w) 



R(u,M) 



(21) 



Jdc{u) - Jdc(-U) 

for each disorder configuration in a ring with M sites with Jdc(w) = nMj and jdc from 
Eq. (jTUl) . The distribution of this rectification parameter is, on symmetry reasons, an even 
function of u, hence (R(u,M)) = 0. In the case of self-averaging, the variance {R(u, M) 2 ) 
should decrease as ~ 1/M for M — » 00. As shown in Fig. [3l this behavior is nicely confirmed 
by taking disorder averages of R 2 (u, M). 
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FIG. 3: Variance (R 2 (u, M)) of the distribution of the rectification parameter R(u,M) for the 
ring system in dependence of the system size M at two fixed values of the bias u. The R(u, M) 
were calculated from Eqs. (|12l2ip for a box distribution of site energies with A e = 6 and disorder 
averages were performed over 10 3 — 10 5 realizations. 

IV. CONDUCTION IN OPEN CHANNELS 

So far we have considered ring systems with periodic boundary conditions. In many 
situations the coupling of the system to particle reservoirs is of importance, as in molecular 
wires, ion channels through membranes, and thin-film electrolytes in contact with non- 
blocking electrodes. In these systems details of the contact with the reservoir can play 
a decisive role for the transport behavior, so that a specific treatment is needed for the 
particular system under consideration. 

On the other hand, if one is interested in generic features of the particle transport, one 
can adopt a coarse-grained description, where only a few external parameters enter, as, for 
example, the thermodynamic driving force of a reservoir to bring the system into equilibrium 
with itself. Based on such coarse-grained description we will in the following characterize 
a reservoir by its chemical potential (amounting to a "site energy level" relative to the site 
energies of the system), and an energy barrier for exchanging particles between the system 
and the reservoir. 

To be specific, we consider a one-dimensional channel consisting of M sites, which is 
coupled to sites k = and k = M + 1, belonging to two reservoirs with chemical potentials 
/i° = eo an d /Ur = caz+i, respectively. Particles are injected or ejected from the two reservoir 
sites with rates that fulfill the condition of detailed balance with respect to the grand- 
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canonical ensembles associated with /i£ and /i R . As for the ring system, the site energies e k 
and the barrier energies Uk,k+i, k — 0, . . . , M, determine the jump rates in the absence of 
the external bias u, see Sec. [II] (U 0i i and Um,m+i specify the energy barriers for exchange 
of particles with the left and right reservoir, respectively). In the presence of a spatially 
uniform bias u, the potential drop along the channel leads to the site energies 

E k = e k -ku, (22) 

and the electrochemical potentials 

^ = E Q = nl and /x R = E M+1 = fi° R - (M + l)u , (23) 

if we locate the point of zero external potential at the left end of the channel. Note that for 
k = and k = M, Eqs. ( TTH314I) define the jump rates for entering and leaving the system, in 
agreement with detailed balance with respect to the grand-canonical ensembles associated 
with /x° = e and /i R = €m+i- 

In the open channel the particle number is a random variable and it is not possible to 
consider a single-particle approach from the beginning. The rate equations for the local con- 
centrations Pi = (rii) follow from a Fermi lattice gas model, where the occupation numbers 
rii at each site can have only two values = (vacant site) or rii = 1 (occupied site), and 
the set {rii] specifies the microstate in the channel. The average (...) has to be taken with 
respect to the probability distribution of the microstates at time t, whose time evolution 
follows a master equation. Based on the master equation the derivation of the currents jij+i 
in the equations of motions ([5]) is straightforward (for a systematic approach, including also 



models with particle-particle interactions going beyond site exclusion, see Ref. 1281 ). The 
result is 

j i)i+1 = r+(n i (l-n m ))-rr +1 (n m (l-n i )) ) i = 1, . . . , M - 1 . (24) 
For the boundary currents specifying the exchange of particles with the reservoirs one obtains 



3o,i = Tq (1 - pi) - T 1 p 1 , 

jm,m+i = ^mPm — r M+1 (i — Pm) . 



(25a) 
(25b) 
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In a mean-field approximation, (n$ni+i) — (ni)(n i+ i) = PiPi+i, the currents in Eq. (I24|) 
can be expressed as 



= r+pi(i - - r i+1 p i+ i(i - pi) 



z = 



1,...,M-1. 



(26) 



In contrast to the ring system, the pk are no longer normalized, but the mean number density 
p of particles is, for fixed energy disorder, controlled by the electrochemical potentials /il 
and /ir. Accordingly, the currents jk,k+i i n Eqs. ( |26f25l) are particle currents (rather than 
probability currents) along the bonds between sites k and k + 1. 

The nonlinear dependence on the pt leads, for non- vanishing bias u > 0, to interesting 
phase transitions of the mean particle concentration with respect to variations of /zl and 
/xr, even in systems without energetic disorder.^ Based on exact solutions of the nonlinear 
mean-field rate equations, one can show that these phase diagrams are correctly predicted 
by the mean-field approximation.— The fact that phase transitions can occur also in the 
dilute limit is sometimes disregarded. For example, it has not been considered in treatments 
of incoherent hopping transport of electrons along DNA molecules. 

A thorough study of the nonlinear Eq. (1261) in the presence of energetic disorder goes 
beyond the scope of this work. In the special case of pure barrier disorder (all q = 0) and 
a current driven solely by a chemical potential difference A/i = /i° — /^r (bulk bias u = 0), 
one has Tf = T^ +1 , and the nonlinear terms cx PiPi+i in Eq. ( 1261) cancel. Accordingly, an 
analytical solution of Eqs. (I25|26l) can be obtained for the stationary state following the 
procedure discussed in the following Sec. IIV Al The result for the corresponding dc-current 
reads 



Note that due to the physical meaning of the jk,k+i discussed above, the (total) current J^ c 
appears in Eq. (1271) . 

In the further treatment we will focus on situations where the consideration of the 
one-dimensional geometry is an approximation for a preferred bias direction of a higher- 
dimensional system, i.e. the pi in Eqs. ( T251I261) are mean concentrations (per site) that rep- 
resent averages over a larger number of sites belonging to lines or planes perpendicular to 
the current direction. In this case we can, without worrying about the boundary-induced 




(27) 



(r+ + 17) + (r+ + r- +1 ) + (r + + 17) (r+ + r M+1 ) £ 
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phase transitions in one- dimensional geometries, consider the dilute limit of Eq. (126]) with 
1 -p% ^ 1, 

Ki+1 = Tf Pi - rr +lPi+1 , i = 1, . . . , M - 1 . (28) 

The rate equations for the occupation probabilities p k (t) now have the same form as in 
Eq. (jSJ) for the single-particle transport on the ring, but we have to take into account the 
boundary currents according to Eq. (12"5"|) . Moreover, one should keep in mind that the pi, 
according to the derivation of Eq. (I26I) . should be much smaller than one.— 

In total five external parameters control the transport behavior in our model for the 
open channel: The chemical potentials /x° and /i R , the energy barriers ZTo.i and Um,m+i for 
particle exchange of the system with the reservoirs, and the bias u. In the following, we 
will in most cases consider the (J>r, ^o,i, Um,m+\ to be given and discuss the transport 
behavior with respect to the bias u. 



A. DC current 



To calculate the stationary current under a static bias we iterate Eq. (jSJ) to obtain 

pt = A k _ lP f - J dc B k ^ , (29) 

with 

k 

A k = Hvi = exp[(£i - E k+1 )} , (30) 
i=i 

k k k 

B k = ^2 K rn Yl Vl = ^2 Km eX P( E m+l ~ E k+l ) , (31) 
m=l l=m+l m=l 

where the expression containing the products hold true in general, while the second expres- 
sions are valid for detailed balanced rates. 

Using Eqs. (1251) together with Eq. (|29|) for k = M one obtains a closed equation for J dc 
with solution 

, l-exp[-(/i L -^)-(M+lH 

1 exp[-e Af - Mu - /i°]) exp(e fc - - ku) 

r+ r+ 2-^ r+ 

1 M k=1 1 k 

= LzIEEtM , 02b) 

exp(-Efc - )Ul) 



r + 

fc=0 fc 
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where A/i = A/io — (M + l)u. This result in turn fixes the local concentrations pf in Eq. (129]) 
via Eqs. (I301I3B and pf from Eq. fl25ij) with j ,i = ^de- 
Equation (132bj) may be interpreted in a similar way as the linear response in the ring 
system, cf. Eq. ffT3"l) : The current follows from a driving force 1 — exp(— A/i) and a to- 
tal "conductance" given by elementary "conductances" exp[— (E k — fih)]^t m ser i a l order. 
Equation (132bl) is, however, not a linear response formula, but describes the full nonlinear 
response to the bulk driving force u and the boundary driving force A/io- Note that these 
driving forces do not enter Eq. (132bl) in the single combination A/i = A/i — (M + l)u, since 
and are controlled independently by fJ[ and /z^, respectively. 

Due to the factors oc exp(— ku) in Eq. (I32al) . only jump rates (r^ /fe ), k = 1, 2, . . ., from 
sites close to the left (right) boundary give a significant contribution for positive (negative) 
bias u. This means that for u ^ 0, Jdc is governed by jump rates belonging to sites in a 
region of size oc 1/u close to either boundary. As a consequence, already pure barrier disorder 
(with all €k = 0) leads to rectification effects in the open channel, in marked contrast to the 
behavior in the ring system. 

It may be surprising at first sight that the dominant contribution to the current comes 
from regions close to either boundary (for similar phenomena expected in connection with 
electron transport though molecular bridges, see Ref. |3l|). The effect can be understood 
when considering, without generality, u > 0, and a single large barrier > Uq in an 

otherwise ordered system with smaller barriers Uk,k+i — Uq for k ^ I (and all = 0). Let 
us first look at the density profile in the region of sites left [k < 1} and right [k > (I + 1)] 
of the large barrier. For the current jij+i > across the large barrier Uij + i to equal all 
other currents jk,k+\i the concentrations p& in the right region have to be much smaller than 
Pj, while the local concentrations in the left region must decrease smoothly with increasing 
distance from the large barrier (smaller k) . Hence the density profile in the stationary state 
has a maximum at site j with a smooth decay to the left and a sharp fall to the right of the 
large barrier Ujj + i. This is demonstrated in Fig. HI where we show the solution pf 1 for a large 
barrier close to the left boundary (solid line) and close to the right boundary (dashed line). 
As a consequence, when the large barrier is closer to the left boundary, the density at the 
boundary site k = 1 becomes larger, leading to a smaller current Jd c = jo,i- More generally 
speaking, we can say that for u > (u < 0) the energy landscape close to the left (right) 
boundary controls the density at the boundary site k = 1 (k = M) and thus the current 
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^dc = jo,i (Jdc = 3m,m+i)- We note that the dominance of the boundary regions will no 
longer apply when considering the transport with site exclusion in strictly one-dimensional 
topologies (ASEPs or TASEPs). 



3x10" 
2x10" 
1x10" 

°0 5 10 15 20 25 30 

k 

FIG. 4: Stationary density profiles in an open channel with M = 30 sites, a constant bias u = 1, 
and a single large barrier Uij+i = 5 close to the left (/ = 4, solid line) and close to the right 
boundary (I = 26, dashed line); the other barriers are set to one, U^k+i = 1 for k ^ I (including 
the boundary barriers for exchange of particles with the reservoirs with ^ R = \ft L = —10), and all 
efc = 0, k = 0, . . . , M + 1. As a consequence of the density profile, the current J = 1.8 • 10 -5 for 
the large barrier at site / = 4 is smaller than the the current J = 4.7 • 10~ 5 for the larger barrier 
at site I = 26. In the latter case J has practically the same value as in the corresponding ordered 
system (all U k ,k+i = 1)- 

To illustrate typical behaviors of the current, we calculate Jdc as a function of the driving 
forces for only barrier disorder (all e k = 0) and for only site energy disorder (all Uk^+i = 0), 
using the box distributions introduced in Sec. [Til Figures show results for the disorder 
averaged current (a) as a function of u for //° = /i^ = —10, and (b) as a function of A/io for 
pP = —10 (/^l r = — 10 ± A/io) and u — 0. Similar as in the ring system, the current-voltage 
curves in Fig [5^ have a more convex shape in the presence of site energy disorder for small 
u. One may ask if the current j™ 8 in the ring system [Eq. ( fl2l) ] and the current in 
the open channel [Eq. (132j) ] can be connected by simply taking account the mean number 
N = Yl!k=iPk of particles in the channel, i.e. if = NJ™ g . However, the fact that regions 
close to either boundary govern the value of J^, already shows that such mapping cannot 
be correct. Indeed, based on the analytical results (j!2|32p obtained for the ring system and 
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open channel, one can show that such a relation does not hold true. Numerical solutions also 
show that the relation does not provide a reasonable approximation (see also the discussion 
in Sec. ED- 




FIG. 5: Current J^ c (a) as a function of the bias u at fixed fjS = ^ = —10, and (b) as a function 
of the chemical potential difference A/j,° for vanishing bias u = 0. Averages have been performed 
over 100 realizations of the disorder, for a box distribution of energy barriers with Ajj = 5, and a 
box distribution of site energies with A e = 6. Solid lines refer to the exponential rates and dashed 
lines to the Glauber rates. In the case of barrier disorder and u = 0, the exponential and Glauber 
jump rates are the same, and hence the corresponding currents agree in part (b). 



B. Thermodynamic limit and rectification 

The dominance of the boundary regions implies that the thermodynamic limit has to 
taken is such a way that for u > the left boundary has to be fixed and the right boundary 
goes to infinity, while for u < one should consider the reversed situation (fixed right 
boundary and left boundary going to infinity). We focus on the case u > here (with 
obvious analogous treatment for the case u < 0). For M — >• oo, Eq. (I32a[) then becomes 

E exp(e fc - ku - /i L ) >^ exp(E k - /i L ) 

k=0 k k=0 k 

One can proof that for point-symmetric energy landscapes (e& = eM+i-fe, £4,fc+i = 
Uu+i-k M-k) the current is antisymmetric with respect to a reversal of the driving forces, i.e. 
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Jdc{—u, fih — > /iR,/iR — > fih) = — J(u, A*r) (the reference point of zero external potential 
has to be shifted from the left to right boundary also). 

Moreover, as mentioned above, rectification effects occur already for pure barrier disorder 
and do not become smaller for increasing M. Accordingly, the width of the distribution of 
the rectification parameter defined in Eq. (T2T1) should saturate to a finite value for M — > oo. 
This is confirmed in Fig.[6j where for pure energy disorder, {R 2 (u, M)) is shown as a function 
of M for two fixed values of u and /i° = /x R = —10. It would be interesting to check this 
theoretical prediction in experiments, e.g. in thin film electrolytes contacted to non-blocking 
electrodes. Systematic measurements in dependence of the system size (film thickness) would 
allow one to distinguish between a possible finite size effect and the effects induced by the 
open boundaries. 



0.12 F 




FIG. 6: Variance (R 2 (u, M)) of the distribution of the rectification parameter R(u,M) for the 
open channel in dependence of the system size M at two fixed values of the bias u. The R(u, M) 
were calculated from Eqs. (|32l2ip for a box distribution of site energies with A e = 6 and disorder 
averages were performed over 10 3 — 10 5 realizations. 



V. TIME-DEPENDENT NONLINEAR RESPONSE 

In this section we discuss the time- dependent nonlinear response to a sinusoidal electric 
field E(t) = Eq sin(u;t) with large amplitude E , corresponding to a bias u(t) = -u sin(cG>t) 
with amplitude uq = qE a/k-QT > 1. To this end we solve the rate equations ([5]) supple- 
mented by periodic boundary conditions for the ring and Eqs. (125!) for the open channel. 
After a transient time interval the stationary regime is reached, where we determine the total 
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current J s t(t) averaged over many periods. Fourier decomposition of this stationary current 
yields the complex first order and higher harmonics J n {uo) = J' n {ui) + iJ"(u), n = 1, 2, . . . 

In the high and low frequency limits the current J s t{t) (and hence the harmonics J n (uj)) 
can be calculated analytically. For u — > oo and barrier disorder, the mean local densities 
Pi(t) in the stationary state become independent of position and time, i.e. Pi(t) = p,— and 
one can show that for each realization 

[/+(«(*))-/-(«(*))], (34) 

where f±(u) are the factors modifying the transitions due to the external driving [see 
discussion before Eq. f ll5p ]. Upon averaging over the disorder (or due to self-averaging), 
Y^h=i ex P( — Uk,k+i)/M can be replaced by the ensemble average (exp(— C/i )2 ). 
For uj — > 0, one can take the quasistatic limit, 

J s t{t) = Jdc (u(t)) (35) 

with Jdc(-) from Eq. ffl2|) for the ring system and Eq. fl32|) for the open channel. For 
exploring the intermediate frequency behavior we have to rely on our numerical solution of 
the underlying rate equations. 

In the following we will concentrate on barrier disorder, implying that harmonics of even 
order vanish in the ring due to the absence of rectification (see the discussion in Sec. IIII Bj) . 
In the open channel, by contrast, rectification effect are present and the harmonics of even 
order are nonzero. However, these harmonics of even order are much smaller than the 
harmonics of odd order, and therefore will not be shown here. For the discussion of the 
harmonics of odd order we focus on the real parts J^(^). 

Figure [7] shows the harmonics J[{oj) and J%{oj) in the case of the exponential jump rates 
for the barrier disorder with AU = 2 and bias amplitude Uq = 1 (for the channel we have 
set /i° = = —1 and boundary barriers U 0t i = Um,m+i = 2.2. The results were averaged 
over 5 realization of the disorder. The circles mark the results for the ring system and the 
squares for the open channel. 

In the ring system, the first harmonics J[(u)) shows the typical behavior known for a hop- 
ping system in the linear response limit: In a high frequency regime, J[(uj) shows a plateau, 
and then, upon lowering the frequency, it decreases monotonously within a dispersive regime 
until approaching the low-frequency regime, where J[(u>) again becomes independent of u. 
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FIG. 7: First order and third order harmonics of the current in the ring and open channel for 
Afj = 2, M = 2000, and the exponential jump rates (/xl = [J-r = — 1 and Uo,i = Um,m+i = 2 for 
the open channel). The results have been averaged over the same sets of 5 realizations of barrier 
disorder, and the single-particle results for the ring are matched to the mean particle number in 
the open channel. The dotted lines mark the limiting behavior for high and low frequency (see 
text). 

The third order harmonics J^(u>) in the ring also shows a plateau at high and low frequencies, 
and passes through a minimum in the dispersive regime. The plateau values in the limits of 
high and low frequencies follow from Eq. (134j) and Eq. (135j) . respectively, and are marked by 
dotted lines in the figure. With respect to the imaginary parts J"(oj) and J£{uj), we found 
peaks appearing in the dispersive regimes in Fig. [7J 

In the open channel the harmonics follow those in the ring system for higher frequencies. 
This can be understood from the fact that at higher frequencies the dynamics in the interior 
of the channel is dominant ("bulk behavior"). At lower frequencies, however, the coupling 
to the reservoirs leads to significant changes in the mean particle number. As a consequence, 
an additional dispersive regime^ is seen at low frequencies, until the limit corresponding to 
Eq. (j35p is reached. Note in particular that J'^ioS) changes its sign when approaching the 
low-frequency limit. 

Let us finally note that we have obtained an analogous overall behavior of the harmonics 
in the case of site energy disorder with the notable difference that no change of sign in J^(lj) 
was observed. 
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VI. SUMMARY AND CONCLUSIONS 



The problem of one-dimensional hopping transport has gained renewed interest, in par- 
ticular in connection with biophysical applications and electron transport through molecular 
wires. We have discussed in this work the situation for non-interacting particles with a focus 
on disorder effects (or regular variations of site and barrier energies) on the current response 
to an external bias. For both the periodic ring system and the open channel analytical 
results were derived for the stationary current in response to static external driving forces, 
without making specific assumptions on the form of the jump rates. Representative results 
were shown for spatially uncorrelated energy landscapes, characterized by box distributions 
either in the barrier or site energies. 

It was further shown that in the ring system rectification effects become smaller for 
increasing system size. In the thermodynamic limit of infinite system size, the current Jdc(w) 
becomes anti-symmetric with respect to the bias u and its expansion in powers of u can 
exhibit non-analyticities of the form |ii| 2n+1 w, n — 0, 1, ... In the open channel rectification 
does not vanish in the thermodynamic limit due to the fact that the current is dominated 
by the variations of the energy landscape close to either system boundary dependent on the 
bias direction. It would be interesting to check this rectification effect in experiments, as, 
for example, in measurement of ionic currents in electrolytes in contact with non-blocking 
electrodes. 

Numerical solutions of the underlying rate equations were obtained for a sinusoidal exter- 
nal driving and results were presented for the first and higher harmonics of the current. For 
intermediate and high frequencies the harmonics in the open channel were shown to equal 
those in the ring, if the particle concentration is adapted properly. In the low-frequency 
regime the harmonics can be derived from the quasistatic limit. This implies that the low- 
frequency limit is different in the open channel from that in the ring. The origin of this 
difference can be attributed to changes in the mean particle number in the open system, 
which are not present in the ring model. 

The results presented here provide a basis for further investigations of interacting parti- 
cles. As discussed in Sec. IIV} in truly one-dimensional geometries already hard-core inter- 
actions can change the general characteristics of the transport behavior due to boundary 
induced phase transitions of the mean particle concentration. Influences of disorder effects 
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on these phase transitions have been discussed in various works (see e.g. Refs. I34|j35lj36l ). but 
a thorough general treatment for arbitrary disorder has not been provided yet. Only a few 
studies have been performed for longer range particle-particle interactions. An example is 
the treatment of nearest-neighbor repulsions in TASEPs on the basis of specific rules for the 
transition rates.— ^ This can give rise to more complex phase diagrams compared to the 
case of hard-core interactions. A more complete exploration of the effects of disorder and 
particle-particle interactions, as required to get a more detailed description of real systems, 
still remains an open challenge. 
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